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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 


ROCKET-MODEL INVESTIGATION OF LONGITUDINAL STABILITY AND 
DRAG CHARACTERISTICS OF AN AIRPLANE CONFIGURATION 
HAVING A 60° DELTA WING AND A HIGH 
UNSWEPT HORIZONTAL TAIL 


By Robert F. Peck and Jesse L. Mitchell 
SUMMARY 


A rocket-propelled model of an airplane configuration having a | 
60° delta wing with an NACA 65A003 airfoil section апа a high unswept `. 
horizontal tail has been flown by the Langley Pilotless Aircraft Research 
Division. The results obtained at Mach numbers between 0.9 and 1.37 
indicate little variation in lift-curve slope and in static and dynamic 
stability with change in Mach number and no abrupt trim changes. Sta- 
bility parameters and lift-curye slope show effects of varying lift, 
however, which might be explained by the horizontal tail being in a posi- ` 
tion where the downwash variation with angle of attack increased with ` . 
lift. At & Mach number of approximately 0.90, the model was pitched to 
angles of attack above 20°. This probably resulted from reduced’ stabil- 
ity, at the higher lift coefficients, stemming from the high tail loca- 
tion. The exposed wing lift-curve slope showed little variation with: 
lift and agrees satisfactorily at supersonic speeds with values calculated 
by an approximate linearized theory. 


` INTRODUCTION ` 


The longitudinal stability, lift, and drag characteristics of air- 
plane configurations at transonic and supersonic speeds are being inves- 
tigated by the NACA using rocket-propelled models in free flight. The 
effects of variations in wing geometry on the longitudinal aerodynamic 
characteristics are being investigated in a correlated program in which 
various wings are test-flown on a basic fuselage-empennage arrangement 
(refs. 1 to 5). The data from these models are obtained from telemetered 
records of the response of the models to & square-wave variation of I 
horizontal-tail incidence. 
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The purpose of the present paper is to present the results obtained 
from the flight of one of the general research test vehicles which had a 
60° delta wing with an NACA 65A003 airfoil section. Data presented 
include lift, drag, pitching moment, and damping in pitch of the complete 
configuration in addition to exposed wing lift. The Mach number range 6 


covered by the test is from 0.9 to 1.37 and Réynolds number from 8.4 x 10 


to 13.8 x 106. 


The tests меге made at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 

CN normal-force coefficient, ous 
CN normal-force coefficient of exposed wing based on total wing 

M area, Ne/qS τ 
Се chord-force coefficient, — ES 7 
Ст, lift coefficient, Cy cosa ~ Cc sina 
Ст drag coefficient, Со cos a + Cy sin а 
an /& normal accelerometer reading _ 
81/8 longitudinal accelerometer reading 
Ne corrected wing balance reading, lb Е 
g acceleration due to gravity, ft/sec? - _ 
α angle of attack, deg i 
W weight of model, lb 
S wing area (including area inclosed within fuselage) sq ft 
Se exposed wing area, sq ft 5 mE 
q dynamic pressure, lb/sq ft 
Cm pitching-moment coefficient k: PR Жш 
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e angle of pitch, radians 

e angular acceleration in pitch, radians/sec? 

c wing mean aerodynamic chord, ft 

T moment of inerti& in pitch, slug-ft^ 

7 distance between nose and center-of-gravity normal 


accelerometers, ft 


б horizontal-tail deflection, deg 
M Mach number 
ү velocity, ft/sec 
R Reynolds number based on mean aerodynamic chord 
P period of asispa, sec 
11/9 time to damp to one-half amplitude, sec 
t time 
Subscripts: 
ELI 
do c 
q = $t 5 
e exposed 
WE wing fuselage 
Cg center of gravity 


The symbols a, q, and & used as subscripts indicate the deriva- 
tive of the quantity with respect to the subscript. 
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MODELS AND APPARATUS 


Model NUN = ae: I 


А sketch of the configuration used in this investigation is shown 
in figure 1. A 60° delta wing with ап МАСА 65A003 airfoil section was 
mounted on a basic fusehage-empennage arrangement similar to that used 
in references 1 to 5. Füselage ordinates are given in reference 3. 

The lower vertical fin was installed on the model to minimize roll or 
уам motions which might perniciously influence the longitudinal maneuver 
in & manner noted in reference 5. 


The wing was constructed of solid steel; the horizontal tail, of 
solid duralumin; the fuselage, of magnesium and steel; and the vertical 
tails, of wood with duralumin skin. 


The incidence of the horizontal tail was varied in an approximate 
square-wave pattern by means of a hydraulic pulsing system. The tail 
settings used were -1. 16° and -5.439 (with respect to the wing). 


The model weight was 142.3 pounds, the center of gravity was at 
20.6 percent of the mean aerodynamic chord, and the moment of inertia 
in pitch was 11.56 slug-ft2. = 


Instrumentation τ _ 


The model was equipped with a lO-channel telemeter transmitting 
continuous measurements of normal acceleration near the center of gravity 
апа at a point about 3 feet ahead of the center of gravity, longitudinal 
and transverse accelerations, angle of attack, horizontal-tail position, 
two measurements of exposed-wing normal force, total pressure, and a 
reference static pressure. = 


A vane -type angle-of-attack indicator which had a range of 15° to 
~10° relative to the sting was mounted on a bent sting 50 as to provide 
angle-of-attack instrument range up to approximately +21° 


Wing normal force was measured by means of a beam-type wing balance 
equipped with an inductance-type pickup. Gaps at the wing-fuselage 
juncture were sealed with fabric. A strain gage was attached to the wing 
balance primarily to gain experience in telemetering strain-gage measure- 
ments. The wing-normal-force data presented herein were obtained from 
the inductance pickup record. 


The total-pressure tube was mounted on & strut below the fuselage 
as shown in figure 1. А static-pressure orifice was located on the top 
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of the fuselage 0.70 body diameter behind the forward station of the 
cylindrical portion of the fuselage. 


SCR 58} radar was used to determine the flight path of the model 
as a function of time. Atmospheric conditions at the time of the flight 
were determined from radiosonde data. 


Mach number of the model was obtained from the telemetered total 
pressure measurements and the static pressures obtained from the SCR 584 
radar and radiosonde data. 


TESTS AND ANALYSIS 


The model was launched at a 60° elevation angle from a mobile 
launching platform (fig. 2). A 6-inch-diameter solid-fuel ABL Deacon 
rocket motor was used to boost the model to & Mach number of 1.39. The 
model, which had no sustainer rocket, separated from the booster at 
rocket burnout by virtue of its lower drag-weight ratio. The data were 
obtained during the decelerating portion of the flight as the model 
responded in pitch to the square-wave variation of horizontal-tail 
incidence in a manner indicated by the portion of telemeter record shown 
in figure 3. There were no appreciable transverse accelerations during 
the portion of flight over which data were obtained. 


The short-period-pitch oscillations resulting from the deflection 
of the horizontal tail were analyzed by the methods of appendix A, refer- 
ence 1, to determine the trim, longitudinal stability, lift, and drag 
characteristics of the model. Tn addition, tbe data from the two normal 
accelerometers were used to obtain total pitching-moment data as follows: 


The pitching acceleration 6 is given by 
ё = Е (а/е) ове т (an/8) се| 
which is proportional to total-pitching-moment coefficient 


Т oe 
 ntotal ^ qSc š 


The Reynolds number of the test varied with Mach number as indicated 
in figure h}. 
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CORRECTIONS AND ACCURACY 


Corrections —-- 


The wing normal forces obtained from the wing balance included not 
only aerodynamic normal forces but also inertial forces exerted by the 
wing and the moving parts of the wing balance. The total normal forces 
read by the balance were corrected for these inertial forces through use 
of normal-acceleration data and the weights of the contributing components. 
The aerodynamic normal force was converted to coefficient form and the 
exposed wing lift coefficient was assumed given by Che = CNe сов a. 


The angle of attack at the center of gravity was obtained from the 
angle-of-attack measurements made at the nose of the model by the шее 
outlined in reference 6. 


By use of the Cmq + Сш data obtained from this test and the 
assumption that the value of Cma was primarily that contributed by 


the horizontal tail, the relatively small pitching moments due to damping 
were estimated and subtracted from the total pitching moment, ав obtained 
by the method briefly described in the foregoing dde to obtain 

pitching moment due to angle of attack. 


Accuracy 


Possible systematic errors in the absolute values of directly 
measured quantities are proportional to the total range of the instrumen- 
tation involved. Possible errors obtained by assuming an accuracy of 
approximately tl percent of the total instrument range are presented in 
the following table: | = 


Experience has shown the accuracy of variations in the various 
quantities far exceeds the accuracy of the absolute values indicated by 
the foregoing table. This is verified by the relatively small amount 
of scatter of individual points in the basic data plots shown herein. 
The Mach number is believed accurate within 20.01 throughout the test. 
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Ihe angle-of-attack indicator may be subject to a further possible 
error in absolute angle of attack of 10.50 due to small asymmetries in 
the vane which may cause it to float at small angles relative to the 
air stream. 


In this particular investigation, at Mach numbers below 1.15, the 
angle-of~attack measurements were subject to further error as evidenced 
by recorded flat peaks on the angle-of-attack oscillations. These flat 
peaks and a lag in angle-of-attack readings are believed due to excessive 
friction in the indicator which was greatly reduced at the higher Mach 
numbers by shaking of the nose section. This is illustrated in figure 3 
where the portions of telemeter record taken just before and after the 
nose shaking ceased are shown. 


RESULTS AND DISCUSSION 


Trim 


The trim lift and angle-of-attack characteristics for the two tail 
deflections used are shown in figure 5. The solid portions of the lines 
indicate regions where trim was actually measured as the mean line of 
the pitch oscillation. No abrupt trim changes are indicated in the Mach 
number range covered. 


Lift 


Plots of variation of lift (both total and exposed wing) with angle 
of attack during the first one and one-half cycles of the pitch oscilla- 
tions are presented in figure 6. The Mach number change occurring during 
the one and one-half cycles shown was the order of 0.04 for the high lift 
and about 0.02 for the low lift oscillations. The Mach numbers quoted 
in the figure were the average during the portions of oscillations 
presented. 


The difference between Ст, obtained under conditions of a 
increasing with time and those with a decreasing with time results, 
it is believed, primarily from an a lag effect caused by excessive 
friction in the angle-of-attack indicator as discussed briefly in the 
section entitled "Accuracy." The slopes do not appear to be appreciably 
affected by the sign of da/dt except at the peaks of the oscillations. 


The lift-curve plots for а Mach number of 0.90 indicate а lift break 
at an angle of attack of about 14.59. This lift break was apparently 
associated with the wing since it was indicated by the wing balance 
record as well as by the normal accelerometers. А similar but less 
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abrupt break in lift curve is indicated by wind-tunnel data on a delta 
wing of aspect ratio 2 (ref. 7). Differences in indicated abruptness 
of the Lift break are probably due to the comparatively low number of 
wind-tunnel test points in the lift-break range. 


The model pitched beyond the lift break to angles of attack exceeding 
the range of the angle-of-attack indicator and, during the time it was 
at these high angles, the Mach number dropped to values too low to be 
measured accurately by the instrumentation used. Ав the angle-of-attack 
instrument limits were approached, a combined pitch-lateral maneuver 
occurred as indicated by relatively high values of lateral acceleration; 
therefore no data were obtained below a Mach number of 0.90 where the 
lift break and pitch-up occurred. Ж 


Lift-curve slopes were measured about the trim value of Cy, апа a. 


Variations of lift-curve slope with Mach number as obtained for the total 
model and for the exposed wing are given in figures ауа and 7(b), 
respectively. 


Lift-curve slope of the total configuration (fig. T(a)) at the high 
Lifts is indicated to be somewhat lower than that at the lower lifts. 
It is believed that this primarily resulted from the horizontal tail 
being in a position such that the downwash variation with angle of attack 
de/da was greater ав the lift increased. This type variation of de/da 
with Cy, is indicated by data in references 8 and 9. The exposed wing- 
lift-curve slope (fig. 7(Ъ)) indicated only a small amount of nonlinearity 
over the lift range covered, which supports the foregoing statement con- 
cerning the downwash variation. 


Results of previous work such as those given in reference 10 have 
shown that exposed wing-lift-curve slope may be estimated in many cases 
by multiplying the lift-curve slope of the wing fuselage by the ratio 
of exposed to total wing area. The lift-curve slope of the wing fuselage 
used in the present test was obtained by subtracting the lift-curve slope 
of the tail from the lift-curve slope of the total configuration at the 
low lifts. The lift-curve slope of the tail was obtained from data of 


reference 2 and downwash values based on data of references 9 and li. It ` 


can be seen from figure 7(b) that, in this particular instance, using 

the area ratio to approximate wing-lift-curve slope results in values 
appreciably lower than those measured. This indicates that this estimate 
is not applicable to configurations with comparatively low values of Se/S. 


A closer estimate of the exposed wing-lift-curve slope may be obtained 
for supersonic speeds from the approximate linearized theory of 
reference l2. 


There was no definite indication of wing buffeting except near a 
Mach number of 0.90 at lift coefficients at and above the lift break 


— 
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where & relatively mild but irregular shaking was experienced. The 
vibration of the nose which occurred &bove Mach numbers of about 1.15 
(fig. 3(a)) is not attributed to wing buffeting. The frequency of this 
shaking corresponds to & body-nose bending frequency determined by shake 
tests made before the flight апа the occurrence of such а vibration has 
been noted in the same Mach number region on another model with the same 
fuselage-empennage but & different wing plan form (unpublished data). 


Drag 


The basic drag data are presented in the form of polars in figure 8. 
The effect of varying Mach number on drag coefficient in the drag-break 
region is evident in the variation of points in the polar for an average 
Mach number of 0.99. Variation of drag coefficient with Mach number at 
several constant lift coefficients and the two tail settings is given in 
figure 9. The increase in drag coefficient as the Mach number increased 
from 1.2 to 1.35 is contrary to the trends indicated by data of refer- 
ence 5. The reasons for this increase are not known but may be due to 
effects of interference or to the nose shaking shown in figure 3 which 
was evidenced &t Mach numbers above 1.15. А large portion of the change 
in drag coefficient resulting from & change in tail incidence is the 
change in the streamwise component of the tail normal-force coefficient 
with change in tail incidence. 


and with Mach 


acy 2 OT Seis 
number are presented in figure 10. The data indicate that, at low lift 
coefficients (5 = -1.16), the resultant force vector due to angle of 
attack was tilted forward. However, at the higher lift coefficients the 
vector was essentially normal to the wing as indicated by coincidence of 


Variations of the induced drag factors 


acp = Ж]! 
and —————. This effect is in agreement with data presented in 
ас 2 57. 3Chg, 


reference 13. 


Static Stability 


The variation of pitching-moment coefficient with lift coefficient 
is shown in figure 11. The values of Cy presented in this figure were 


obtained from two normal accelerometers by methods discussed in foregoing 
sections of this paper. The variations of period of the pitch oscilla- 
tions and the stability parameters Cm, and ἄσπι/ᾶσι, with Mach number 


are shown in figure 12. Values of dCm/dCy, shown in figure 12(c) were 
obtained from the data of figure 11 and also from Cm, in figure 12(b) 


ΜΗ --- 


10 МАСА RM L52KOLla 


divided by Cz, shown in figure T. Agreement of ἆσπι/άσι, values deter- 
mined from these two methods is good. ΜΙ. - 


The fairly smooth variation of static stability with Mach number was 
similar to variations obtained from somewhat similar delta-wing configura- 
tions reported in references 5 and 13. The stability decreased with 
increase in lift coefficient. This decrease is evident in figure 11 at 
Mach numbers of 0.99 and 0.90. At a Mach number of 0.90 the data were 
obtained almost to a Lift coefficient of 1.0 and indicate a very severe 
loss of stability and a pitch-up at lift coefficients above 0.75 which, 
is the point where the sharp break in lift occurred. 


It is believed that the primary factor contributing to the decrease 
in stability with increase in lift, at lift coefficients below the lift 
break, was an increase in de/da at the teil with an increase in lift. 
This uns associated with the high tail position and was mentioned in the 
section entitled "Lift" of this report as a possible cause of the reduc- 
tion in total-lift-curve slope with increase in lift coefficient. 


Damping in Pitch 


Figure 13 presents the information obtained on the damping charac- 
teristics of the pitch oscillations resulting from the abrupt control 
movements. Figure 13(а) shows the amplitude ratio measured from the 
trim line as a function of time. The data have been faired by a straight 
line on the semilog plots. The time for the oscillations to damp to one- 
half amplitude as determined from the faired curves is^shown in fig- 
ure 13(b). These times to damp, along with the lift-curve slopes from 
figure T, have been used to determine the rotary-damping parameter 


Cmg + σπα (fig. 13(с)). Ап 1псгеаве 1п Сша + Cy. with increase in - = 


lift coefficient is indicated. An increase in de/da “with increase in 
lift coefficient mentioned in the "Lift" and "Static Stability" sections 
of this report could account for this increase in Cmq + σπις through 


its effect in increasing (πὰ. 


Comparisons 


Comparisons are made in figure 14 between data from the present test 
and data from references 2, 5, and 13. Reference 2 contains data on the 
Wingless fuselage empennage used on the model of the present test. The 
longitudinal cbaracteristics of an airplane configuration with the same 
fuselage-empennage group without lower vertical tail and with & modified 
delta wing are presented in reference 5 and the longitudinal character- 
istics of a tailless configuration with a 600 delta wing having an ` 
NACA ©5 (96 4006.5 airfoil section are presented in reference 13. For 
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purposes of comparison, the drag of one vertical tail, obtained from 
reference 2, was added to the drag of the model of reference 5. 


In general, the parameters compared show the same variation with 
Mach number, and differences in level may be explained by differences 
in configurations. The model of reference 13 had no horizontal tail and 
therefore would be expected to differ from the models of the present 
test and reference 5 as indicated. The small differences in,lift-curve 
slope, aerodynamic center, and induced drag between the models of the 
present test and reference 5 are compatible with the slightly higher 
aspect ratio of the model of reference 5. However, the differences in 
supersonic minimum drag of these two models are contrary to what might 
be expected. These differences are in general of the same order as the 
accuracy of the drag data; however, the drag of the present model may 
have been increased by unfavorable interference of the lower vertical 
tail or by nose shaking (at the higher Mach numbers) as discussed briefly 
in the "Drag" section of the present report. 


CONCLUSIONS 


A rocket-propelled model has been flown to determine the Longitudinal 
stability and drag characteristics of an airplane configuration with a 
60° delta wing with an NACA 65A003 section and a high unswept horizontal 
tail. The data indicate the following: 


l. The variations of lift-curve slope, static stability, and pitch- 
damping parameter (mq + Cms) with Mach number were smooth and relatively 


small and no abrupt trim changes were encountered. 


2. Total-lift-curve slope and static stability decreased and the 
pitch-damping parameter (Cm, + Cm. ) increased with increasing lift coef- 


ficient, all of which could be accounted for by the horizontal tail being 
in a position where there was an increase in rate of change of downwash 
with angie of attack as the lift coefficient increased. 


3. The exposed wing-lift-curve slope showed only a small variation 
with change in lift coefficient, agrees satisfactorily at supersonic 
speeds with values calculated by an approximate linearized theory, but 
is appreciably higher than values estimated by multiplying lift-curve 
slope of the wing fuselage by the ratio of the exposed to total wing 
area. 


h. At & Mach number of approximately 0.90, the reduced stability at 
the higher lift coefficients, stemming primarily from & high tail loca- 
tion, resulted in the model pitching up to angles of attack above 20°. 
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An abrupt break in lift was indicated both by the total lift and wing 
lift measurements at a total lift. coefficient of approximately 0.75. 


5. In general, data obtained at low lift show the same effects of 
Mach number variation as comparable data from other delta-wing tests. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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Figure 1.- General arrangement of the model. All dimensions are in inches. 
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Figure 2.- Photograph of model and booster on launcher. Ne d 
ЄС | κα. 
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(a) During nose shaking. 
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(b) After shaking ceased. 


Figure 3.- Portion of telemeter record taken during апа after nose shaking. 
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(а) Trim angle of attack. 
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(р) Trim lift. 


Figure 5.- Trim lift and angle-of-attack characteristics. 
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(b) Exposed wing. = 


Figure 6.- Variation of lift coefficient with angle of attack. 
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Figure Τ.- Lift-curve slope. 
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Figure 8.- Drag polars. 
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Figure 9.- Variation of drag coefficient with Mach number апд. 
i | lift coefficient. 
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Figure 10.- Induced drag factor. ` 
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Figure 11.- Variation of pitching-moment coefficient with Lift coefficient. 


αν 


z 
г 
Ξ 
М 
B 
3 


Gc 


26 | CC EET NACA RM L52KO0la 


P,sec —— TJ 
«8 1.4 
М 
(а) Period of pitch oscillation. 
-.04 
na -.02 
O 
.8 
М 
(b) Longitudinal stability parameter Св. 
л Cma/CLa 
оп 40/00] from figure 11 
ас _/ac, 
.8 .9 1.0 τα 1.2 1.5 1.4 
М 
(ο) Longitudinal stability parameter aC, (С. 
Figure 12.- 


Static longitudinal characteristics with center of gravity 
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(a) Variation of amplitude ratio with time. 
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Figure 13.- Damping-in-pitch characteristics. 
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— Present test 


Test of reference 5 
plus a vertical tail 
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(a) Lift-curve slope. 


(b) Minimum drag coefficient. 


Aariynamia center, 
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(d) Induced drag factor. 


Figure 1%.- Comparisons of low-lift data (5 = -1.16°) with data of othe 


tests. 
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